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Abstract 

Exclusive and inclusive (z/^/x - ), (u e ,e~) cross-sections and /^"-capture rates 
are calculated for 12 C and 16 O using the consistent random phase approxima- 
tion (RPA) and pairing model. After a pairing correction is introduced to the 
RPA results the flux-averaged theoretical (z^,/i~), (z/ e ,e~) cross-sections and 
/i~-capture rates in 12 C are in good agreement with experiment. In particular 
when one takes into account the experimental error bars, the recently mea- 
sured range of values for the (zv„,/i~) cross-section is in agreement with the 
present theoretical results. Predictions of fi~) and (v e ,e~) cross-sections 

in 16 O are also presented. 
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Investigations of neutrino-nucleus interactions are presently in the forefront of particle 
and nuclear astrophysics research. The study of basic properties of the neutrino such as 
its mass, possibility of decay or oscillations, is probed using the interaction of neutrinos 
with nuclei. The neutrino interaction with matter is very weak and therefore neutrinos that 
are emitted from the inner parts of a star are "messengers" that bring information about 
the processes occurring inside the stars Q. These extraterrestrial neutrinos are detected in 
reactions involving nuclear targets. It is of considerable importance to provide a theoretical 
framework that takes into account the various aspects of nuclear structure and is able to 
describe well these reactions . Attempts of this kind were made in the past [[J and in 
the last few years such attempts have been revived 0,|J in view of the availability of new 
experimental results in this field of neutrino-nucleus interaction. In particular we will refer 
here to the KARMEN results and the results obtained at the Liquid Scintillator Neutrino 
Detector (LSND) in Los Alamos. The 12 C(i/ M , /i _ )X measurements at the LSND were made 
with the same experimental setup as the one used in the neutrino oscillation experiment - 
the results of which were recently published 0. A study of the 12 C( t / M ,/i~)X reactions will 
shed indirectly some light on the question of validity of the neutrino oscillation experiments. 
For example a strong disagreement between theory and the experimental results for the 
12 C(z/ At , /i~)X could probably lead to some skepticism concerning the experimental setup of 
the neutrino experiments in general. In fact, the recently published theoretical results 0] 
for the inclusive 12 C(z/ M , /i - )X are in disagreement, by a factor of two, with the LSND result 

In the present paper we examine a range of neutrino -nucleus cross-sections as well as 
capture rates with a special emphasis on the 12 C case. Our calculations are performed first 
in the Hartree-Fock (HF) approximation and then in the consistent Hartree-Fock- Random 
Phase Approximation (HF-RPA). The consistency refers here to the fact that the HF mean 
field and the particle-hole (p-h) interaction result from the same effective nucleon-nucleon, 
two-body force. Such RPA calculation preserves the energy-weighted sum rule of the strength 
distribution of one-body operators @ and therefore it is a favorable scheme when it comes to 
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calculate the distribution of total strength. One is led to expect that, for closed shell nuclei 
inclusive cross-sections of processes governed by one-body transition operators should be well 
reproduced in an RPA-type calculation. Indeed, the RPA had much success in providing 
an adequate description of giant resonances |8||| and of a variety of inclusive processes 
in nuclei with good closed shells. An example relevant to our subject of neutrino-nucleus 
interactions is the inclusive //"-capture on nuclei. The HF-RPA was able to reproduce |§ the 
inclusive capture rates in a number of closed shell nuclei. It was found that the collectivity 
and RPA ground state correlations are very important in reaching good agreement with 
experiment. This suggests that other inclusive neutrino processes such as for example the 
(i>^,fi~) reaction (which is the inverse to /i-capture) will be well described by the HF-RPA. 
It is known however that the ground state of 12 C is not a good closed shell. Admixtures 
of the (pi/2) 2 G°3/2) 2 configuration into the (^3/2) 4 configuration are large and one should 
expect substantial corrections to the matrix elements obtained in the RPA. As we will see 
these corrections play a very important role when one calculates the (1/1, l~) cross-sections, 
in particular the exclusive ones to the ground state of the daughter nucleus. The inclusive 
cross-sections to the excited states are less affected but, in order to achieve agreement with 
experiment one must nevertheless include these corrections also for the excited states. 



The numerical applications are performed with the Skyrme forces SGII and SIII fllQllTl 



The interaction SGII was adjusted so as to give the correct value of the Landau parameter 



Fq JXTj] in the spin-isospin particle-hole channel which should be of particular importance in 
charged current neutrino reactions on nuclei. First, the Hartree-Fock equations are solved 
in coordinate space to obtain the self-consistent mean field. This mean field determines 
the single-particle spectrum. For the present problem it is not necessary to treat exactly 
the single-particle continuum since we are not studying specific (exclusive) channels where 
a nucleon would be emitted following the {p^yT) or (u e ,e~) reaction. Therefore, it is 
convenient to discretize the single-particle spectrum by diagonalizing the Hartree-Fock mean 
field on a harmonic oscillator basis. The reference Hartree-Fock state \HF) corresponds 
to the target nucleus, e.g., 12 C. Let us denote by (a, &,...) the proton (neutron) 



occupied states and by I, J,... (A,B,...) the proton (neutron) unoccupied states. The 
proton (neutron) creation and annihilation operators are respectively pf and Pi (nf and 
m). In reactions of the (z/^/z - ) or (u e ,e~) type the final states |A) belong to the ATz= -1 
daughter nucleus (e.g., 12 N) and they can be described by the charge-exchange RPA [pT2|JT3[ 
model: 

I A) = (E^ff Pfna + E^A Ptn A )\0) , (1) 

I, a i,A 

where |0) is the correlated RPA ground state. The and are solutions of the charge- 



exchange RPA equations |12| , ^L3| . For a one-body charge-exchange operator of the general 
form: 

O = E Oa0pin fi , (2) 
the transition amplitude (A|F|0) can be expressed simply as: 

(A|O|0) = 52xi?0 Ia -Y. Y iA*O iA . (3) 

la iA 

In the case of a parent nucleus with zero angular momentum in the ground state the 
cross-section is [@|§: 

a = ^-cos\Qc)Y.ViEiT{Z,E l ) t d(cos(6))M X0 , (4) 

Z7T A J-l 

where G and 9c are the Fermi constant and the Cabibbo angle, pi and Ei are the momenta 
and energies of outgoing leptons (muon or electron), 9 is the angle between the momenta 
of the lepton and the incoming neutrino. The factor T accounts for the effects of the final 
state interaction (FSI) of the outgoing lepton with the daughter nucleus of charge Z [0,0. 
For the case of {y ei e~) reactions the mass of the outgoing lepton is small and the effect is 
not so important. The effect of the FSI for the negatively charged muon is more significant, 
increasing the cross-section approximately by 15 — 20%. In Eq.([|), the sum goes over the 
available nuclear excitations, denoted by A. The nuclear structure effects are incorporated 
into A^ACb the bilinear combination of the nuclear matrix elements between the ground state 
|0) and the excited states |A) of the daughter nucleus. These are given by [§]: 



Mxo = M F |(A|F|0)| 2 + M GT |(A|GT|0>| 2 + M' GT A (5) 

The coefficients M* are obtained by the Foldy-Wouthuysen transformation of the weak 
Hamiltonian where the terms up to third order in the momentum transfer q/M are kept 
(M is the nucleon mass) ||. The first matrix element squared of Eq. fl5|) is: 

|<A|F|0>| 2 = 47rX;|(A I J||t_j J (?r)i r j||0> | 2 (6) 

Here, <f is the momentum transfer, q — \q\, a and £_ refer to the nucleon spin Pauli matrices 
and isospin-lowering operator, respectively, ||...|| stands for the standard definition of the 
reduced matrix elements, ji are the spherical Bessel functions and Yj are the spherical 
harmonics. The remaining combinations of the matrix elements are: 

|<A|GT|0>| 2 = 47r£ |(A, J\\t_ 3l (qr) [Y t x ||0)| 2 , (7) 
i, j 

A = 4tt f-V /2 y (-l)^')/ 2+J ((2/ + l)(2l' + I)) 1 ' 2 

(I I' 2 ) ( 1 1 2) 

\(\,J\\t-.ji(qr)\Yix<?\ J \\0) 
{ 1 2 J { V I J J 

x (A,J|M/'(^)[^'X^ljl|6)*. (8) 

Here, [ x ]j means the coupling to the total angular momentum J. In our calculation 
we take into account states with J < 3 with positive and negative parity 

We perform the integration over angle 9 in Eq.(|j) with the step size A# = 27r/30. The 
single-particle matrix elements of the operators O are calculated using the Hartree-Fock 
wave functions in steps of Ar=0.1 fm in the radial coordinate. 

The 12 C nucleus is not described well by a closed subshell and configuration mixing 
is important in the ground state. Besides the RPA correlations, one of the important 
correlations is introduced by the pairing force. Here, we estimate the effect of pairing on the 
exclusive and inclusive neutrino cross-sections. In the expressions (f§-[8[) for the cross-section, 
two types of single-particle matrix elements enter: a) those which do not contain the spin 
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operator a and b) those which contain a. The two kinds have different symmetry properties 
under time reversal. Consequently, the corresponding expressions for these single-particle 
operators in terms of quasiparticles are ||14j| : 

(i*\Oa\k*) = (uiV k - UkVi)(i\Oa\k) , 

{i*\0\k*) = (u^ k + u k Vi)(i\0\k) , (9) 

where the asterisks mean that the pairing has been taken into account, O is an operator 
that depends on r only, and the u and v are the coefficients of the Bogoliubov-Valatin 
transformation. In the case of 12 C, the u and v factors differ significantly from their values 
without pairing only for the p%/2 and pi/2 levels. This results in reduction factors: 

C = ( U Pl/2 V P3/2 ~ V Pl/2 U P3/2) ' (10) 

which multiply the single-particle matrix elements of the operators containing a in the 
expressions for the RPA matrix elements (A|O|0). The evaluation of pairing effects in the 
cross-sections is done by introducing the factors u^k ± u k v , in the matrix elements obtained 
from the RPA solutions without pairing. 

To obtain the cross-sections that allow for comparison with the experimental data one 
has to fold the energy-dependent cross-section of Eq.(|J) with a corresponding neutrino flux 
f(E) 

(a) f = J dEa(E)f(E) , (11) 
where f(E) is a properly normalized neutrino flux from an available neutrino source, 

hE) = m^m • (12) 

and f{E) is the initial (unnormalized) flux from the source. Here, the value of E depends 
on the neutrino source used in each experiment. It is taken to be zero for the case of the 



electron neutrino l|15|, while for the case of the muon neutrino source Eq = E t \ ir pi, where 



E thr is the threshold energy for the A Z{y ili , \i ) reaction. We calculate the cross-sections with 



the neutrino energy steps AE = lMeV and AE = 5MeV for the electron and muon neutrino 
cases, respectively. The spectra for the muon and electron neutrinos were taken from Refs. 
P]T5|j. The endpoint for the electron neutrino flux is 52.7MeV, thus the cross-section for 
the A Z(z/ e ,e~) reaction is sensitive to the low-energy excitations in the daughter nucleus, 
mainly the giant resonances. For the case of the muon neutrino, the flux is cut at E= 260 
MeV. 

The results of the calculations are shown in Table 1. We note first that the results 
obtained do not differ very much for the two interactions used and the cross-sections or 
capture rates agree typically within 10% for the SGII and SIII interactions. In order to 
assess the influence of collective effects and of the RPA ground state correlations we have 
calculated each cross-section or capture rate first in the Hartree-Fock approximation and 
then in the HF-RPA. The comparison shows that in the RPA the /i~-capture rates are 
reduced by 30% and 50% in 12 C and 16 0, respectively. This is in agreement with the results 
in Ref. 0. Note that in the present calculation of the //"-capture rates the SU4 assumption 
is not used and the vector, axial vector and induced pseudo-scalar contributions are directly 
calculated. From Table 1 we see that the RPA effects are also very important in the (z/ e , e~) 
process. In 12 C and 16 O the flux averaged (z/ e , e~) cross-section is reduced by nearly a factor 
two when one goes from HF to HF-RPA. 

One of the manifestations of the fact that 12 C is not a closed shell is the need for a 
large suppression factor in order to reproduce the experimental exclusive cross-sections to 
the ground states of the A=12 nuclei @||. In Ref. [§] a reduction factor C 2 =3.88 - 4.13 
was introduced in the computation of exclusive processes, however these authors have not 
considered the influence of configuration mixing in the 12 C ground state on the inclusive 
cross-sections to the excited states. In the present work, by including the pairing correlations 
we are able to treat the effects of such configuration mixing in both exclusive and inclusive 
cross-sections. The pairing calculations are done employing Eq. (|) and a t> p i/2=0.65, a 
choice in accordance with the value of the gap parameter A ~3-4 MeV |T4J. The pairing 



correction reduces the flux-averaged exclusive cross-sections and /i -capture rates to the 



ground state of the daughter nucleus by factors 4-7 compared to the RPA results, depending 
on the type of reaction one calculates. The inclusive cross-sections to the excited states 
(denoted as < a* >) are affected less by the pairing correction but still the effect is sizable, 
reducing the < a* > cross-sections by 25%. The reduction due to pairing of the total in- 
clusive processes is about a factor of 5 for the (u e , e~) cross-section(because it is dominated 
by the transition to the J 77 = 1 + ground state) and about 33% reduction for the 
cross-section with respect to the RPA result. Note that the RPA and pairing correction 
decrease the cross-sections calculated in the HF approximation by a factor of 6-7 for (z/ e , e~) 
and a factor of 2 for the (v^, fi~) and /i"-capture processes. These reductions are essential 
in achieving agreement with experiment. We emphasize that the agreement is achieved in 
all the quantities calculated by using the same value of the parameter v p i/ 2 - We did not 
attempt to find the best value for this parameter. In parenthesis of column 4 we show the 
results obtained when we use t> p i/2=0.60 . We should remark that the results in Table 1 
corresponding to interaction SGII are obtained with a slightly changed spin-orbit parameter 
(compared to the original SGII force) so as to reproduce in the RPA the experimental thresh- 
old energy. Our results for the (z/^,/1 - ) inclusive cross-section of 13.5xl0 _40 cm 2 (for SGII) 
and 14.5xl0 _40 cm 2 (for SIII) should be compared with the recently revised experimental 
value of (11.2±0.3±1.8)xl(T 40 cm 2 from the LSND experiment @. 



In 16 O the nucleons form a good closed shell and pairing effects are not important. Except 
for the ^-capture there are no experimental neutrino data. Our HF-RPA predictions in 
16 for the flux-averaged inclusive (u^, ff~) cross-sections are 27. 8xl0 -40 cm 2 (SGII) and 
27.1xl(r 40 cm 2 (SIII) whereas for (i/ e , e") they are 17.2xl0~ 42 cm 2 (SGII) and 16.9xl0" 42 cm 2 
(SIII). The calculated HF-RPA /i"-capture rates are 1.05xl0 5 s -1 (SGII) and 1.0Qxl0 5 s _1 
(SIII) to be compared with the experimental value of 0.98xl0 5 s _1 ||17|| . 

In summary, we have shown that the use of the consistent HF-RPA scheme and the 
introduction of pairing corrections can successfully reproduce the experimental yU~-capture 
rates and neutrino cross-sections in 12 C. Although the present approach is not as elaborate 
as an extended shell-model calculation or a quasi-particle RPA it contains the basic features 
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that emerge from such more extensive methods. 
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TABLES 

TABLE I. Flux-averaged cross-sections and //"-capture rates in 12 C. a exc stands for exclusive 
cross-sections to the 12 N ground state, a* is the inclusive cross-section to excited states. The 
capture rates oj g , s , and totot correspond to the partial rate to the ground state and to the total rate, 
respectively. The (z^,/i~) cross-sections are in 10~ cm 2 , (u e ,e~) cross-sections are in 10 -42 cm 2 , 
capture rates are in 10 4 s _1 . For each quantity, values calculated with SGII and SIII are in the 
upper and lower row, respectively. The results in brackets correspond to the choice 1^1/2=0.60 . 



Channel 


HF 


HF 
+pair. 


HF- 
RPA 


HF-RPA 
+pair. 


Exp. 




4.28 
4 70 


0.61 
fi7 


3.35 
3 80 


0.39 (0.64) 
50 (0 77) 


0.66±0.1±0.1 |U 


a* 


22.7 
24.1 


16.1 
18.0 


17.7 
18.6 


13.1 (13.7) 
14.0 (14.4) 






27.0 
28.8 


17.7 
18.7 


21.1 
22.4 


13.5 (14.3) 
14.5 (15.2) 


11.2±0.3±1.8 @ 


@exc 


78.1 
100.4 


10.0 
14.2 


54.8 
68.2 


7.1 (11.4) 
10.1 (16.0) 


10.5±1.±1. |18,1£ 
8.2±0.65±0.75 |1 


)| 


a* 


8.6 
14.0 


8.8 
10.7 


8.3 
8.1 


5.6 (6.2) 
6.4 (6.7) 


6.4±1.45±1.4 [|C 

3.6±2.7 || 
5.7±0.6±0.6 [21 
8.6±1.2±1.5 [15 


1 
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0",V. 



90.6 
114.4 



18.8 
24.9 



63.2 
76.3 



12.9 (17.6) 
16.5 (22.7) 



14.1±2.3 IH 

i6.8±i.7 m 



// -capt. 



ff-s. 



3.61 
3.70 

8.0 
8.4 



0.51 
0.53 

3.87 
4.10 



3.24 
3.48 

6.87 
7.22 



0.41 (0.67) 
0.45 (0.73) 

3.09 (3.48) 
3.23 (3.64) 



0.60±0.04 |22|,£23| 



3.7±0.1 
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